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We show that in Gauss-Bonnet gravity with negative Gauss-Bonnet coefficient and 
without a cosmological constant, one can explain the acceleration of the expanding 
Universe. We first introduce a solution of the Gauss-Bonnet gravity with negative 
Gauss-Bonnet coefficient and no cosmological constant term in an empty (n + 1)- 
dimensional bulk. This solution can generate a de Sitter spacetime with curvature 
n(n + l)/{(n — 2)(n — 3)|a|}. We show that an (n— 1) -dimensional brane embedded 
in this bulk can have an expanding feature with acceleration. We also considered 
a 4-dimensional brane world in a 5-dimensional empty space with zero cosmologi- 
cal constant and obtain the modified Friedmann equations. The solution of these 
modified equations in matter-dominated era presents an expanding Universe with 
negative deceleration and positive jerk which is consistent with the recent cosmo- 
logical data. We also find that for this solution, the ll n"th derivative of the scale 
factor with respect to time can be expressed only in terms of Hubble and deceleration 
parameters. 



I. INTRODUCTION 

The idea of brane cosmology has attracted a great deal of attention recently. It may 
provide us a new solution for the so-called hierarchy problem or the source of dark energy 
and dark matter . This theory is also consistent with string theory, which suggests that 
matter and gauge interaction (described by an open string) may be localized on a brane, 
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embedded into a higher dimensional spacetime, while the field represented by closed strings, 
in particular gravity, propagate in the whole of spacetime. In the brane world scenario, we 
live in a 4-dimensional (4-D) hypersurface embedded in a higher dimensional bulk spacetime. 

Among many brane models, those proposed by Randall and Sundrum which are motivated 
by superstring/M-theory [^] are the most applicable ones. In their first model (RS I) Q], 
they proposed a mechanism to solve the hierarchy problem with two branes, while in their 
second model (RS II) 5], they consider a single brane with a positive tension, where 4-D 
Newtonian gravity is recovered at low energies even if the extra dimension is not compact. 

Brane world cosmology underscores the need to consider gravity in higher dimensions. 
In this context one may use another consistent theory of gravity in any dimension with a 
more general action. This action may be written, for example, through the use of string 
theory. The effect of string theory on classical gravitational physics is usually investigated 
by means of a low energy effective action which describes gravity at the classical level 
This effective action consists of the Einstein-Hilbert action plus curvature-squared terms and 
higher powers as well, and in general gives rise to fourth order field equations and bring in 
ghosts. However, if the effective action contains the higher powers of curvature in particular 
combinations, then only second order field equations are produced and consequently no 
ghosts arise J7j. The effective action obtained by this argument is precisely of the form 
proposed by Lovelock [8]. The appearance of higher derivative gravitational terms can be 
seen also in the renormalization of quantum field theory in curved spacetime 

These facts provide a strong motivation for investigating how the usual features of brane 
world cosmology are modified by more general gravitational theories such as Lovelock gravity. 
In this context, many authors are extensively studying the brane world cosmology in Gauss- 
Bonnet gravity jioj . In these analysis, one needs the solutions of Gauss-Bonnet gravity in 
the bulk spacetime. Because of the nonlinearity of the field equations, it is very difficult to 
find out nontrivial exact analytical solutions of Einstein's equation with hi ghe r curvature 
terms. However, a few exact solutions of this theory have been obtained 

On the other hand, it seems established that at the present epoch the Universe expands 
with acceleration instead of deceleration along the scheme of standard Friedmann model. 

n 

This follows directly from the observation of high red-shift supernova [13( and indirectly 
: rom the measurement of angular fluctuations of cosmic microwave background fluctuations 
1^ | . The latter shows that the total mass/energy density of the Universe is very close to 
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critical value one (Q = 1), while the observations of the Universe in large scale structure 
indicate that normal gravitating (visible and invisible) matter can contribute only 30% into 
the total one. Thus, one concludes that the remaining 70% is some mysterious agent that 
creates the cosmological acceleration. The simplest suggestion is that the source of this 
acceleration is the vacuum energy (cosmological constant). However, it meets the two well 
known cosmological and coincidence problems [3]. A second approach is the consideration 
of cosmological components as a dynamical term. This scheme, usually called quintessence, 
can be achieved by adding a scalar field into Einstein gravity. Several forms of potential 
achieve quintessence prescriptions but non of them seems to be directly related to some 
fundamental quantum field theory. A third approach is to obtain a comprehensive model 
derived from some effective theory of quantum gravity which through an inflationary period 
results into the today accelerated expansion of the Universe. From this point of view, 
gra vitational theories including higher order curvature terms naturally come into the game 
|l6l |. However, some of these modified gravitational theory have their own problems [17| . 

Here we want to explain accelerating Universe by use of the Gauss-Bonnet gravity without 
a cosmological constant term or a scalar field or any other kind of dark energy. Thus, we 
restrict ourself to the second and third terms of Lovelock gravity. The second term is the 
Einstein-Hilbert term, while the third term is known as the Gauss-Bonnet term. From a 
geometric point of view, the combination of the Gauss-Bonnet terms constitutes, for five- 
dimensional spacetimes, the most general Lagrangian producing second order field equations, 
as in the four- dimensional gravity where the Einstein-Hilbert action is the most general 

n 

Lagrangian producing second order field equations |18| . 

The outline of our paper is as follows. We give a brief review of the field equations in Sec. 
ITT1 In Sec. II I II we introduce a solution of the Gauss-Bonnet gravity with negative Gauss- 
Bonnet coefficient and no cosmological constant in an empty (n + l)-dimensional bulk, and 
show that this solution can generate a de Sitter spacetime. In Sec. HXlwe first introduce the 
modified Friedmann equations, and then obtain the solutions of these equation in matter- 
dominated era, which can discuss today accelerated expansion of the Universe. We finish 
our paper with some concluding remarks. 
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II. FIELD EQUATIONS IN GAUSS-BONNET GRAVITY 

The most fundamental assumption in standard general relativity is the requirement that 
the field equations be generally covariant and contain at most second order derivative of 
the metric. Based on this principle, the most general Lagrangian in arbitrary dimensions 
is the Lovelock Lagrangian. The Lagrangian of the Lovelock theory, which is the sum of 
dimensionally extended Euler densities, may be written as 

1 [n/2] 

i=i 

where q is an arbitrary constant, [x] denotes the integer part of x, and Li is the Euler 
density of a 2i-dimensional manifold, 

r _ I _n\-ixaibi...a,ibisr> cidi <i-> Cjdj /<->\ 
V L ) u c 1 d 1 ...c i d i '^a 1 b 1 •••• /v a i b i • l z I 

In Eq. (J2J TZabcd is the Riemann tensor and is the totally antisymmetric product 

of i Kronecker deltas, normalized to values and ±1. Here the first term C\L\ = TZ, is 
just the Einstein-Hilbert, and the second term c 2 C 2 = aiJZ^sTZ^' 5 - AK^RT + TZ 2 ) 
gives us the Gauss-Bonnet term. Of course, one may add a constant term to the above 
Lagrangian, playing the role of cosmological constant term. But, as we mentioned before 
this creates its own problems and therefore we don't disturb ourself with it. The Lagrangian 
of Gauss-Bonnet gravity with zero cosmological constant is 

£-g = \{R + a(n, uy5 n^ 5 - m^-RT + n 2 )}, (3) 

where TZ, TZ^ upa) and 1Z^ U are the Ricci scalar and Riemann and Ricci tensors of the space- 
time, and a is the Gauss-Bonnet coefficient with dimension (length) 2 . Here we restrict 
ourselves to the case a < 0. 

Let us suppose, as would be in the co-dimension one brane world scenario, that the 
(n + l)-dimensional bulk, A4 is split into two parts by a hypersurface S, whose two sides 
will be denoted by S±. Their normal unit vector, n M , will be taken to point away from the 
surface and into the adjacent space. The gravitational action for the spacetime (A4, g pu ) 
with boundary (£, h a b) can be written as 



Ig = 7^2 [ dx n+1 y^C G dx n V^h [K + 2a(j- 2G ab K ah )] , (4) 

^ JM K JY.± 
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where K^ v is the extrinsic curvature of the hypersurface S with induced metric h ab , G ab is 
the Einstein tensor of the metric h ab and J is the trace of 

J ab = l - (2KK ac K\ + K cd K cd K ah - 2K ac K cd K db - K 2 K ab ) . (5) 

The first term in Eq. © is the Einstein-Gauss-Bonnet terms, and the second term is the 
boundary term which is chosen such that the variational principle is well-defined This 
term is, in fact, the generalization of Gibbons-Hawking boundary term in Einstein gravity. 
If we also include a matter contribution to the action 



1 mat 



- / dx n+l ^C^ - [ dx n V^hC%\ (6) 

J M JT, 

then the variation of the total action I = I G + J mat over the metric tensor g^ u gives the 
equation of gravitational fields in the bulk as 

-2nn^ + ak^ii x u + 4n pa n^ ua - 2n; aX n vpaX \ = n% u , (7) 

where T^ v is the energy momentum of the bulk defined by T^ v = 25C ( ^ 4 ^ /5g fJ-u — g^Lm 
Assuming the Z2-symmetry for the brane, then the variation of total action with respect to 
h ab gives Q 

K ab - Kh ab + 2a(3J ab - Jh ab + 2P acbd K cd ) = -—T ab , (8) 

where T ab = 25CS ] /5g ab - g ab C^ ] is the energy momentum in the hypersurface and P ac bd is 
the divergence free part of the Riemann tensor of the metric h ab 

Pacbd = Racbd + 2h a [ c Rd]b + 2h b [dR c ] a + Rh a [ c hd]b- (9) 

Also it is worthwhile to mention that the energy-momentum conservation on the hypersurface 
may be written as 

D b T ab = -2% v h/n^ (10) 



where D b denotes the covariant derivative corresponding to h, 



ab ■ 



III. THE DE SITTER SOLUTIONS IN AN EMPTY BULK 



Here we want to obtain the (n + l)-dimensional solutions of Eq. (J7J) in vacuum, which 
generates a de sitter spacetime. We assume that the (n + l)-dimensional spacetime be 
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isotropic and homogenous, i.e. it has a maximally symmetric n-dimensional space. Since we 
are interested in cosmological solutions, we take a metric of the form 

dr 2 



ds 2 



-dt 2 + f 2 (t) 



1 — kr 2 



+ rW 



(11) 



where dfl 2 is the metric of an (n — l)-sphere and A; is a constant which can take only the 
value 1, and —1. To find the function F(t), one may use any components of Eq. (JJJ). The 
simplest equation is the tt-component of these equations which can be written as 



(f 2 + k) \ a \(n-2)(n-3)(f 2 + k)-f 



0. 



(12) 



where the overdot denotes a derivative with respect to the t coordinate. The solutions of 
Eq. (|12|) can be written as 



/(*) 



2y/AB 



Aexp 



t 



v /(n-2)(n-3)|a| 



+ B k | a | exp 



^(ra-2)(n-3)|a| 

(13) 

where A and B are two arbitrary functions. Also for negative k, fit) = y/—kt is another 
solution. The function fit) given in Eq. (fT3*j) shows that the in + l)-dimensional bulk 
spacetime has an expanding feature. For example, if one choose a brane in the hypersurface 
of r = ro, then the volume of (n — l)-brane is 

2W 2 



V, 



n—l 



r~\t). 



(14) 



r(n/2)' 

It is worthwhile to note that the function f(t) has exponential behavior for the case of flat 
space with k = 0. That is 

t 



f(t)=Aexp 



(15) 



V / (n-2)(n-3)|a| 

Also, one may note that the spacetime is asymptotically de Sitter for an arbitrary value of 
k. In these two cases (k = and t — > oo) the curvature of the spacetime is n[n + l)/{(n — 
2)(n — 3)|a|}, which means that the spacetime is de Sitter. 



IV. MODIFIED FRIEDMANN EQUATIONS IN GAUSS-BONNET GRAVITY 

Now we consider a 5-D bulk spacetime with a single 4-D hypersurface at r = r . The 
most general gravitational equation in five dimension which is symmetric, divergence free 
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and linear in second derivative of the metric is Einstein-Gauss-Bonnet equations. Thus, we 
use these equations in order to consider the 4-D brane world cosmology. The brane world 
(B, h^) is located at the hypersurface r = r , and its induced metric is = g^ u — n^n u , 
where n M = (0, ^/g~^, 0, 0, 0) is the normal unit vector to the brane. Since we assume to have 
a perfect fluid in the brane world at r = ro, the bulk has a spherical symmetric 4-D space. 
Thus, the bulk metric can be written as: 

dr 2 



ds z 



-dt 2 + f 2 (t) 



c 2 (r) 



+ dx + &(%)(# + sin' (pdip 1 



(16) 



where b(x) f° r the closed, flat and open Universe is 



b(x) 



sin x k — +1; 
X k = 0; 

sinh x k = — 1 ; 



(17) 



respectively. The induce metric is 



ds 2 



-dt 2 + a 2 (t) [d X 2 + b 2 ( X )(d(j) 2 + sin 2 0# 2 )] 



(18) 



where a(t) = rof(t). The nonvanishing components of the extrinsic curvature of the hyper- 
surface are 



K xx = c a(t), 

K H = c a(t)b 2 (x), 

= c a(t)b 2 (x)sin 2 (p, 



(19) 



where Co = c(r ). Using Eqs. 
tensor of a perfect fluid 



and (fTTIj) . and the expression of the energy-momentum 



T ab = (p + p)u a u b + ph ab , 
one obtains the following modified Friedmann equations: 

d 4co — n 2 pa 
a 16co|a| ' 

2Ac \a\(k + a 2 ) + lQcl\a\ - 6c a 2 - n 2 pa 3 = 



(20) 

(21) 
(22) 



Equation ()21|) indicates that in matter-dominated era (p = 0), our Universe expands with 
acceleration instead of deceleration along the scheme of standard Friedmann model. To be 



8 



more clear we obtain a class of solutions of Eqs. (}2*T]) and ()22|). First, we go through the 
conservation of energy- momentum. Using Eqs. (fTUj) and ()20|) for the metric (fT6*j) . and the 

fact that the bulk is empty, one obtains 

d( pa 3 ) = -pda 3 . (23) 

Given an equation of state p = p(p), one can use Eq. (}23|) to determine p as a function 
of a. Knowing p as a function of a, one can find a(t) for all time by solving Eq. ()22j) . 
Incidentally, the solution a(t) determined in this way will automatically satisfy Eq. (j21)l . for 
by differentiating Eq. (|2~2~j) with respect to time and using Eq. (J2~3*j). one obtains Eq. (}2~T]) . 



A. Matter dominated era 



Assuming the energy density of Universe is dominated by non-relativistic matter with 
negligible pressure, then p oc a~ 3 , and therefore Eq. ()22|) can be written as 



2Ac \a\(k + a 2 ) - Qc a 2 + lQc 3 \a\ - K 2 p a 3 = °> ( 24 ) 



where po an d cto are the values of p and a at the present. Equation (|2^|) can be solved 
exactly. One can write the solutions of this equation as 



a(t) = Aexp — = + B k exp = . (25) 





where B k = (16cg|a| + 24|a|fcco — k 2 pod^) / (24A) . Since Co is arbitrary, therefore A and 
Bk are two arbitrary constants. The arbitrary constants A and Bk, and the Gauss-Bonnet 
parameter a may be fixed by observational data. In order to give the way of fixing A, Bk and 
a, it is more convenient to work with the Hubble (H = a~ l a), deceleration (q = — iJ _2 a _1 a), 



jerk (j = H~ 3 a~ 1 a a), snap (s = if _4 a _1 a) and ... parameters 2l|. It is remarkable to note 
that for a(t) given in Eq. (|23|) . one encounters only with Hubble and deceleration parameters 
(Ho and qo) at present epoch. Indeed, for a(t) given in Eq. (J25J) d n a/dt n = H n q^ 2 \ and 
therefore the Taylor expansion of a(t) is 

oo 

<t) = a Y,^H n \q \ [n/2 \t-t T, (26) 

n=0 

where [n/2] denotes the integer part of n/2. 
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In order to use the astronomical data for fixing (A, B^, a) or equivalently (do, H Q , q ), 
one may use the relation between the luminosity distance d^ and the redshift z of a luminous 
source. It is easy to show that (see the Appendix): 



di 



z 



1 - 2q - 3q 2 + 



k 



H 2 n 2 



1 

+ 24 



1 - 7q - 2Aql - 15g 3 + 



2fc(l + g ) 

H 2 n 2 
n u Q 



z* + ... 



(27) 



Having enough astronomical data, one can determines (do, H , q ) and therefore (A, Bj., 
a). The solution (|2*5|) indicates that during the matter-dominated era, the Universe expands 
with acceleration instead of deceleration with positive jerk. This feature is consistent with 
the observation of high red-shift supernova or the measurement of angular fluctuations of 
cosmic microwave background fluctuations 14, Also, one may note that for a(t) 
given in Eq. (|2*5j) = ^-^/jaj) -1 and = = —qoo = 1 as t —>■ oo. 



V. CLOSING REMARKS 



In this paper, we added the Gauss-Bonnet term with negative coefficient to the Einstein 
action without a cosmological constant term, and introduced a solution of the field equations 
in an empty (n + l)-dimensional bulk. We found that this solution is de Sitter with constant 
curvature n(n + l)/2|a| for the case of flat space (k = 0) or as t — * oo. We showed that 
an (n — l)-dimensional brane embedded in the bulk can have an expanding feature with 
acceleration instead of deceleration. 

We also considered the 4-dimensional brane world in a 5-dimensional empty space with 
zero cosmological constant. This was done in Gauss-Bonnet gravity which has the most 
general gravitational field equation in five dimension. We obtained the modified Friedmann 
equations and showed that these two equations are not functionally independent if one 
use the conservation of energy. Equation (f2~T]) indicates that our Universe expands with 
acceleration instead of deceleration in matter-dominated era. This was done without need 
to any mysterious fluid with large negative pressure or the cosmological constant. We also 
obtained the solutions of the modified Friedmann equations for the matter-dominated era. 
This solution presents an expanding Universe with positive acceleration and jerk which is 
consistent with the recent cosmological data. We also found that a(t) can be expressed in 
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terms of do, H and go om y> and all the other parameters such as jerk, snap and so on can 
be written in terms of H and go- 

As stated before, the Gauss-Bonnet gravity is the most general gravitational field equation 
in five dimension. In higher dimension one should use more terms for action in Lovelock 
theory. The application of the above method in higher dimension, which needs the use of 
Lovelock gravity with more gravitational terms remains to be carried out in future. 

APPENDIX 



In this appendix, we obtain the relation between the luminosity distance d^ and the redshift 
z. This relation has been derived for Robertson- Walker spacetime in 2l| . Here we obtain 
it for the metric ([18)1 . The physical distance travelled by a photon emitted at time and 
absorbed at the current epoch t in term of z is J^j] 



1 + ^ 
2. 



z + 



2 D 



3 /-, x 5 o 1 . 5 s 
1 + 2*>(1 + Qo) + q% ~ ^0 ~ Y2 qo3 ° ~ 24 



z s + 0(z 4 



while d L can be written as 
d L = a(t ) 



2 r o 



-Kx) 



a(Q a(t - D) v/w a(t - D) L 6 
Recall that for a null geodesic in the spacetime (fTHj) 



{i--x 3 + o( x 5 )}. 



Ao 



r to dt r to cdt r TT . 

~M= — \l + H (t -t) + 



2 + go 



(t - tf 



6(1 + go) + jo u3 

tin 



(to - tf + 0[{to - t) 4 ] 



(28) 



(29) 



°U + \HoD + 
a I 2 



2 + go 



(H Df 



+ 



6(1 + g ) + jo 

24 

Now using Eqs. P5 |) - (j3*U|) . one obtains 



[HoDf + O [(HoD) 4 ] 



dj.iz) 



l h , l 



1 

+ 24 



1 - g - 3gg + j 



kc 2 
H 2 n 2 



2 - 2g - 15g^ - 15gg + 5j + 10g j + s + 



2fcc 2 (l + 3g ) 



_n o 



+0(z 4 



(30) 



(31) 
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